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We theoretically and experimentally investigate the chaotic regime of optical frequency combs generated in 
nonlinear ring microresonators pumped with continuous wave light. We show that the chaotic regime reveals 
itself, in an apparently counter-intuitive way, by a flat top symmetric envelope of the frequency spectrum, when 
observed by means of an optical spectrum analyzer. The comb demodulated on a fast photodiodc produces 
a noisy radio frequency signal with an spectral width significantly exceeding the linear bandwidth of the 
microresonator mode. © 2012 Optical Society of America 
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Oscillators are fundamentally nonlinear systems with 
intrinsically rich dynamics. They exhibit chaotic behav- 
ior, or stable oscillation, based on the choice of their 
operating parameters. In the case of optical oscillators, 
any and all examples of dynamical behavior can be, and 
have been, observed. In particular, onset of oscillation 
in mode locked laser has been known to be associated 
with both hard and soft excitation, determined by the 
dynamics of the system. Soft excitation refers to oscil- 
lation that self starts due to the presence of noise in 
the oscillator circuit. Hard excitation refers to a regime 
whereby an external impulse is required to start the os- 
cillation. In this Letter we investigate both chaotic and 
stable regimes of oscillation in frequency combs gener- 
ated in monolithic ring microresonators made with trans- 
parent materials possessing Kerr nonlinearity pQ. These 
optical combs have been the subject of intense interest 
lately in connection with applications in a variety of ar- 
eas in science and technology [2]. However, a satisfactory 
answer to the question of coherence of these combs, es- 
pecially for those exhibiting a wide frequency spectrum, 
has been elusive. We show that a large number of op- 
tical Kerr combs reported in the literature are in fact 
chaotic, despite their evidently regular and symmetric 
shape. These chaotic combs will naturally display a noisy 
RF spectrum when demodulated on a fast photodiodc, 
as has been reported in various studies. 

The existence of chaotic regime is consistent with the 
intrinsic nature of Kerr frequency combs. What is sur- 
prising is that a comb exhibiting a regular spectral enve- 
lope with a flat top, and with modes symmetrically posi- 
tioned about the pump, can be chaotic. As shown below, 
the observed shape of the comb is in fact the result of 
time-averaged chaos, and is the consequence of the con- 
ventional means for observing these combs, namely the 
use of an optical spectrum analyzer (OS A). 

The observed incoherent frequency comb generation in 
a magnesium fluoride (MgF2) whispering gallery mode 
(WGM) microresonator, and results of numerical simu- 



^10" 
E 

n: 

"S1 1 



I 10" 

u.10 

*10~ 3 
SlO- 4 
S10" 5 



(a) 



CD 



E -10 
20 
l.- -30 
8 -40 
o -50 
5" -60 
3 -70 




20 40 60 80 
Mode number 



100c 



O 



1536 1544 1552 
Optical wavelength, nm 















(c) 




















J 


















f 




\ 
















L 





-12 -9 -6 -3 3 6 9 12 
Normalized spectral frequency, 2nf/y 




-3-2-101234 
Detuning from 34.7 GHz, MHz 



Fig. 1. Optical spectrum of an incoherent Kerr comb 
((a) and (b)) and associated RF signals ((c) and (d)) 
generated by the comb on a fast photodiode. Parts (a) 
and (c) represent the theoretical simulation results, (b) 
and (d) are for the measurement results. To simulate the 
spectrum of the chaotic comb (a) we time averaged the 
instantaneous chaotic comb spectra over t ave = IOO/7. 
This averaging corresponds to the averaging introduced 
by an OSA used to observe the comb spectrum in the 
experiment. 



lation reproducing the experimental data are shown in 
Fig. [1] The envelope has a characteristic flat top dis- 
tinct from envelopes of other combs that also can be 
generated in the system. We found that the RF signal 
produced on a fast photodiode by this comb has a spec- 
tral width that significantly exceeding the linewidth of 
the resonator modes, resulting from the chaotic dynam- 
ics of light confined in the microresonator. The presence 
of chaos is more general than the nonlinear interaction of 
subcombs [3] since it constitutes a fundamental property 
of the nonlinear oscillator described with forced nonlin- 
ear Schrodinger equation [3]. 

The existence of a chaotic comb does not imply that 
the resonator cannot generate a coherent Kerr comb. A 
stable oscillation mode can readily exist in a comb oscil- 
lator and may be observed. This, however, is not trivial 
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since the majority of regular oscillation solutions are as- 
sociated with hard excitation regime [5]. We show that 
an abrupt switching of frequency of the cw pump laser 
facilitates transforming the chaotic comb regime to a reg- 
ular one. 
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Fig. 2. Dependence of the normalized field amplitude 
(A = 651, for the case described in the paper) in the op- 
tically pumped mode on the frequency of the cw pump 
light. Solid lines show stable branches corresponding to 
generation of various Kerr combs (see Fig. [3]), dots cor- 
respond to generation of chaotic combs or combs with 
time-depending envelopes. Plot (a) stands for the soft 
excitation regime [5]. The initial amplitude in this case 
is selected to be IbjKg/j) 1 ^ 2 = 10~ 4 , which approxi- 
mately corresponds to the field in the mode contain- 
ing 1/2 photon. Branches (1-7) are not accessible in 
this case. Plot (b) stands for hard excitation. It is as- 
sumed that ACg/7) 1 / 2 = 10 and b 5 o(gh) 1/2 = 0.3... 3 
and b^gh) 1 ^ 2 = 0.3... 3 to get the solutions. Stable 
solutions are realized for the particular nonzero initial 
condition and detuning selection. The comb can also 
be excited from nearly zero initial excitation of modes, 
if the pump frequency is shifted nonadiabatically, as 
shown in Fig. 01 Constant tuning of the pump frequency 
also allows observing solutions corresponding to stable 
branches [6j[7]. 

We generated a chaotic Kerr comb in a MgF2 WGM 
microresonator. The resonator had 34.7 GHz free spec- 
tral range (resonator radius was approximately equal to 
1 mm), and a loaded Q-factor of about 10 9 . The res- 
onator was produced from a commercially available z- 
cut MgF2 optical window by mechanical polishing. The 
shape of the resonator was that of an oblate spheroid in 
the vicinity of WGM localization. Two evanescent field 
prism couplers were used with the resonator [5]. The 
pump light was injected through one prism and the Kerr 
comb was sampled using the other prism. In this way we 
observed the comb structure excited within the resonator 
and avoided excessive leakage of the pump light to the 



photodiode used to generate RF signal. The resonator 
was pumped with a 1543 nm distributed feedback diode 
laser self-injection locked |9| to an appropriate WGM 
mode. The locking occurred due to the intrinsic Rayleigh 
scattering in the resonators host material. The optical 
power emitted by the laser was approximately 5 mW. 
The input coupling efficiency was approximately 30%, 
and approximately 50% of the light exiting the second 
prism coupler was collected into a multimode fiber and 
sent either to an optical spectrum analyzer, or a fast 
photodiode and an RF spectrum analyzer. 

An example of the observed optical and RF spectra is 
shown in Figs.QJb) andQJd). While the optical spectrum 
is rather regular, the RF spectrum is very wide. It is ap- 
proximately five times wider than the mode bandwidth 
(the bandwidth of the WGMs is approximately 200 kHz). 
This implies that the frequency comb is incoherent. 
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Fig. 3. Simulated spectra of the combs corresponding to 
stable solution branches shown in Fig. [2] 

We performed numerical simulations along the lines 
discussed in [5lH0] . We introduced an interaction Hamil- 
tonian V = —%g(e^) 2 e 2 /2, where g = ?kjgcn 2 / '(Vnjj) is 
the coupling parameter obtained under the assumption 
of complete space overlap of the resonator modes j 1 1 j . 
ojo is the value of the optical frequency of the externally 
pumped mode, c is the speed of light in vacuum, ri2 is 
the cubic nonlinearity of the material, V is the effective 
geometrical volume occupied by the optical modes in the 
resonator, and no is the linear index of refraction of the 
resonator host material. The operator e is given by the 
sum of annihilation operators of the electromagnetic field 
for all the interacting resonator modes that we took into 
consideration e = Y] b-j . 

Equations describing the evolution of the field in the 
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resonator modes are generated using the Hamiltonian 

h = - (7 + + \ [V, Sj] + , (1) 

where Sjo,j is the Kronecker's delta; jo is the number 
of the externally pumped mode, 7 is the half width at 
the half maximum (HWHM) for the optical modes, as- 
sumed to be the same for the all modes involved. The 
external pump is applied to the central mode of the 
group of 101 modes. The pump amplitude is given by 
F = [(2-f 0c P)/{hLo )] 1/2 , where P is the value of the 
power of the pump light. We studied the case of pump 
power given by F(g/j) 1/2 = 10, and anomalous GVD 
set by 2ujq — uj + — w_ = — O.I7. 

The results of simulations is shown in Figs. Ufa) 
andQJc). To obtain the optical spectrum we employed 
the fact that a typical OS A time averages each fre- 
quency component. Therefore, the instantaneous fre- 
quency comb envelope in not the observable in the ex- 
periment, rather it is the comb averaged for times longer 
than 10 ms that is observed. The instantaneous enve- 
lope of the comb looks chaotic and not symmetric, but 
the averaged comb looks regular and similar to the ex- 
perimental observation. The simulated RF signal is also 
similar to the measured RF signal, and is broader com- 
pared with the HWHM of the modes. 

We studied regimes of comb generation numerically 
and found that the chaotic regime is always accessible, 
independent of the initial conditions, Fig. [2] (soft excita- 
tion regime |5j). Generation of regular combs (Fig. |3]), 
on the other hand, depends on the initial conditions 
and cannot be realized if the resonator is initially empty 
(hard excitation regime [S] ) . To achieve steady comb gen- 
eration in an initially empty resonator one needs to nona- 
diabatically change either the power or the detuning of 
the pump light, the two parameters that are accessible 
in the comb oscillator. We demonstrated that an abrupt 
change of the pump frequency result in formation of a 
stable frequency comb Fig. 01 The comb is started from 
chaotic oscillations, but when the laser frequency moves 
to to the position where regular solutions are allowed, 
the oscillation pattern changes. The modification of the 
oscillation regime is accompanied by intricate comb dy- 
namics. This is not the only way for generation of an 
stable comb. It was shown in [6l[7] that scanning the 
pump frequency through the resonator mode also allows 
achieving stable combs. 

To conclude, we have demonstrated, both experimen- 
tally and theoretically, the generation of chaotic Kerr 
combs in a nonlinear microrcsonator. We provide expla- 
nation for why in many experiments the comb spectrum 
is often observed to be symmetric, implying the presence 
of high coherence, while the RF beatnote produced by 
the comb is completely irregular, implying an incoherent 
comb. We have shown that it is possible to tune the comb 
oscillator from the chaotic regime to the regular oscilla- 
tion by a nonadiabatic change of frequency or power of 
the pump laser. 




Time, yt 



Fig. 4. A density plot illustrating time dependence of 
amplitude of comb harmonics. This is an example of exci- 
tation of a regular frequency comb in an empty microres- 
onator via nonadiabatic change of the oscillator parame- 
ters. The comb forms starting at \bj(t = 0)| = 10~ 4 when 
the normalized frequency detuning ((u> — oj )/j) between 
the pump light and the pumped mode abruptly changes 
from —5 to —27 at t = lA/j. A chaotic frequency comb 
is generated if the detuning stays constant at the initial 
level (—5). The abrupt change of the detuning results in 
generation of a regular broad coherent comb at branch 
(1) of Fig. [2j The transient process reveals breather for- 
mation in the resonator. 
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